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We study flavor changing processes T — * B/BX S and J/tp —* D/DX U in the B factories and 
the Tau-Charm factories. In the standard model, these processes are predicted to be unobservable, 
so they serve as a probe of the new physics. We first perform a model independent analysis, then 
examine the predictions of models; such as top color models and MSSM with R-parity violation ; 
for the branching ratios of T — > B/BX S and J/tp — > D/DX U . We find that these branching ratios 
could be as large as 10 -6 and 10 -5 in the presence of new physics. 



I. INTRODUCTION 



The possibility of observing large CP violating asymmetries in the decay of B mesons has motivated the construction 
of high luminosity B factories at several of the world's high energy physics laboratories. These B factories will be 
producing roughly about 10 8 Upsilons. Meanwhile BES has already accumulated 9 x 10 6 J/tp and plans to increase 
the number to 5 x 10 7 in the near future. An interesting question, that we investigate in this paper, is whether the 
large sample of the T and the J/tp can be used to probe flavor changing processes in the decays of T and J/tp. In 
particular we look at the flavor changing processes T — > B/BX S and J/tp — + D/DX U , from the underlying b — *■ s 
and c — ► u quark transitions. For the quarkonium system, these flavor changing processes are expected to be much 
\Q . smaller than in the case of decays of the B or the D meson because of the larger decay widths of the bottonium 
and the charmonium systems which decay via the strong interactions. Indeed the standard model contributions to 
T —>■ B/BX S and J/tp — > D/DX U are tiny. However, new physics may enhance the branching ratios for these 
processes. Whether this enhancement maybe sufficient for these processes to be observable in the next round of 
experiments is the subject of this work. 

Non leptonic decays of heavy quarkonium systems can be more reliably calculated than the non leptonic decays 
of the heavy mesons. A consistent and systematic formalism to handle heavy quarkonium decays is available in 
NRQCD 0] which is missing for the heavy mesons. As in the meson system [|| it is more fruitful to concentrate on 
quasi-inclusive processes like T — > B/BX S and J/tp — > D/DX U because they can be calculated with less theoretical 
uncertainty and have larger branching ratios than the purely exclusive quarkonium non leptonic decays. The branching 
ratios of exclusive flavor changing non leptonic decays of T and J/tp in the standard model have been calculated and 
found to be very small 

We begin with a model independent description of the processes T — > B/BX S and J/tp — > D/DX U . In the standard 
model these decays can proceed through tree and penguin processes. For new physics contribution to these processes 
we concentrate on four quark operators of the type sbbb and uccc. We choose the currents in the four quark operators 
to be scalars and so these operators may arise through the exchange of a heavy scalar for e.g a Higgs or a leptoquark 
in some model of new physics. These four quark operators, at the one loop level, generate effective sb{g,j, Z} and 
uc{g, 7, Z} vertices which would effect the flavor changing decays of the B and the D mesons. The effective vertices 
for an on shell g and 7 vanish and so there is no contribution to b — > sj or c — > wy. We can however put constraints 
on these operators by considering the processes b — ► sl + l~ and c — * ul + l~ . The constrained operators can then be 
used to calculate the branching ratios for T — > B/BX S and J/tp — > D/DX U . 

We then consider some models that may generate the kind of four quark operators described above. A few examples 
of models where these operators can be generated are top color models and MSSM with R parity violation. In some 
cases constraints on the parameters that appear in the prediction for the branching ratios for T — > B/BX S and 
J/tp — > D/DX U are already available. In other cases the parameters are constrained, as in our model independent 
analysis, from the processes b — > sl + l~ and c — > ul + l~ . 
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II. EFFECTIVE HAMILTONIAN. 



In the Standard Model (SM) the amplitudes for hadronic T decays of the type bb — > sb + sb are generated by the 
an effective Hamiltonian ^,|| To the standard model contribution we add higher dimensional four quark operators 
generated by physics beyond the standard model In this paper, we consider the four quark operators with two 
scalar currents. 

L new = §1(1 - J 5 )bb(l + 7 5 )6 + + J 5 )bb(l ~ 7 5 )6 + h.c. (1) 

The four quark operators in L new are the product of two scalar currents. In Eq. (1) A represents the new physics scale 
and R\ and R2 are two free parameters which describe the strength of the contribution of the underlying new physics 
to the effective operators. In our analysis we will only keep dimension six operators suppressed by 1/A 2 and neglect 
all higher dimension operators. The details of the matrix elements for the processes T — > B/BX S and J ftp — > D/DX U 
can be found in 



III. LOW ENERGY CONSTRAINTS AND MODELS 



The lagrangian L new generates, at one loop level, the effective sb"/*, sbg*,sbZ vertices, where 7* and g* indicate an 
off shell photon and a gluon. Similar vertices involving c — > u transitions are generated in the charmonium sector also. 
These vertices, with a 7 and Z , will contribute to b — > sl + l~ and c — * ul + l~ . Note there is no contribution to 6 — * 97. 
The vertex b — > sg* can give rise to the process b — > sqq which will contribute to non-leptonic B decays. We expect 
the constraints from b — > sl + l~ to be better than from non-leptonic B decays because of the theoretical uncertainties 
in calculating non-leptonic decays The additional contribution to the effective Hamiltonian for b — > sl + l~ can be 
written as 

SH b ^ sl+l - = -yj^|| J - x) log (£\ [Risj%lHI + RaS^bah^] ( 2 ) 

which has to be added to the standard model contribution || . Similar results can also be written for the charm sector. 

Now we look at various models that can give rise to L new given in Eq. 1. As a first example we consider a recent 
version of top color models M . In such models the top quark participates in a new strong interaction which is broken 
at some high energy scale A. The strong interaction, though not confining, leads to the formation of a top condensate 
< thtR > resulting in a large dynamical mass for the top quark. The scale A is chosen to be of the order of a 
TeV to avoid naturalness problem which implies that the electroweak symmetry cannot be broken solely by the top 
condensate. In the low energy sector of the theory, scalar bound states are formed that couple strongly to the b quark 

II 

L b = ^=b L (H + iA°)b B + h.c (3) 

where /# ~ 50 GeV is the top pion decay constant. On integrating out the Higgs fields H and A we have an effective 
four fermion operator 
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L ^S = -72^T b LbRbRb L (4) 
J a " l h 

Since the b quark in (4) is in the weak-eigenstate, L e ff in (4) will induce flavor changing neutral current (FCNC) four 
quark operators in Eq. (1) after diagonalizing the quark mass matrix ||, with coefficients, 

R l = 7 f2j,2 \ D Lbb\ D Rbb D* Rbs 

* JTi" L H 

R ' 2 = 7 t 2 ' 1 2 \ D Rbb\ 2 D Lbb D* Lbs (5) 

* Jir m H 
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where Dl and Dr are the mixing matrices in the left and the right handed down sector. In the charm sector similar 
interactions can arise due to the strong couplings of the top quark to top pions. The effective operators generated by 
integrating out the top-pions are similar to Eq. (1) with the replacement of b by c and sbyu. In topcolor II models 
L0|, where there can be strong top-pion couplings of the top with the charm quark, we have 
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~ 72 J, 2 \ U Lc 
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(6) 



In supersymmetric standard models without R parity, the most general superpotential of the MSSM, consistent with 
SU(3) x SU(2) x U(l) gauge symmetry and supersymmetry, can be written as 



(7) 



where Wr is the i?-parity conserving part while W# violates the i?-parity. They are given by 



W$ = A,,,. /.,/. ,/.':. + XUUQjDl + \"UCDW c k + ViLiHi. 



(8) 
(9) 



Here Li(Qi) and Ei(Ui,Di) are the left-handed lepton (quark) doublet and lepton (quark) singlet chiral superfields, 
with i, j, k being generation indices and c denoting a charge conjugate field. H\p. are the chiral superfields representing 
the two Higgs doublets. In the i?-parity violating superpotential above, the A and A' couplings violate lepton-number 
conservation, while the A" couplings violate baryon-number conservation. Ay* is antisymmetric in the first two indices 
and A" jfc is antisymmetric in the last two indices. While it is theoretically possible to have both baryon-number and 
lepton-number violating terms in the lagrangian, the non-observation of proton decay imposes very stringent conditions 
on their simultaneous presence pi] . We, therefore, assume the existence of either L- violating couplings or B- violating 
couplings, but not the coexistence of both. We calculate the effects of both types of couplings. 

In terms of the four-component Dirac notation, the lagrangian involving the A' and A" couplings is given by 



-A' 



i J k 



+ 444 + (4)*(4) c 4 



e L a R u L 



4(4) c 4 + 4(4) c < + u R {di) c d k 



h.c, 
+ h.c. 



(10) 



(11) 



The terms proportional to A are not relevant to our present discussion and will not be considered here. The exchange 
of sneutrinos with the A' coupling will generate L new for T — > BX S with 



Itl — —Z-ii t 

4 m^. 

n ^ A i23 A i33 

4 mi 



(12) 



For the case of J/tp — ► DX U the operators in L new cannot be generated at tree level. 



IV. RESULTS 

The various inputs to our calculations can be found in Q . The standard model contribution to the branching ratio 
is 5.2 x 10~ n from the penguin induced b — > s transition. The process T — » BX S can also have a contribution in the 
standard model from tree level processes. The effective Hamiltonian , suppressing the Dirac structure of the currents, 
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G F 

H w = —=V ub V* s [ai(ub)(su) + a 2 {sb)(uu)] 



(13) 



B + K . We can estimate the branching ratio 



v/2 

where a\ and 02 are the QCD coefficients can generate the process T 
for this process as 

BR[T -» B+K-} « |§-| 2 £fl[T -> 

Using BR[T -> calculated in Ref j| one obtains BR[T -> B+R-] ~ 1.5 x 10~ 14 . For a rough estimate of 

5i?[T -> B+A s ] we can scale BR[T -> B+A"~] by the factor BR[B -> D°X]/BR[B -> Dtt]. The measured value 
of BR[B -> D°A] || includes D° coming from the decay of Z?°* and Z? + *. From the spin phase factors BR[B — > 
L>*A] ~ 3Bi?[B -> DI]. Hence BR[B -> D°X]/BR[B -> L>tt] - 20 leading to Bi?[T -► £+X s ] - 3 x 10~ 13 . 

So far we have not considered i?i and i?2 in L new . In our model independent analysis we vary Ri/A 2 , R2/A 2 one 
at a time and the use the constraint from measurements of b — > se + e~ and b — ► s/j + /i~ . We identify A with the 
masses of the exchange particles which we take to be between 100 — 200 GeV. The allowed values of Ri/A 2 , A2/A 2 
are then used to calculate T — > BX S The constraint from b — ► sZ + Z~ gives 



|fli,a|/A 2 < (6 - 9) x 10-%1/GeV) 2 

Using the upper bounds on |i? 1-2 |/A 2 we find the branching ratio for the process T(15 r ) — > BA S to be between 
(1 — 2) x 10 -6 . Branching ratios of similar order are also obtained for T(25) and T(35). For T(45) the branching 
ratio is smaller by a factor of 100 because of the larger width of T(4S I ) which decays predominantly to two B mesons. 
Turning now to models, we find for the top color model from Eq. (5) we can write 



D 



,Ri 
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A 2 m 2 t \D Lbb \ 



'R>Rbb 



(14) 



We can identify A = ma and use the constraint from b — > se + e~ for a typical value of |i?i|/A 2 ~ 6 x 10 _6 (l/GeU) 2 
. Assuming |-Dl&&| ~ l-D^I w 1, and f% = 50GeV we obtain |Z3flb s | ~ 2m 2 H x 10 -6 . With typical values of 
771 h ~ 100 — 200 GeV we get |Da:& s | ~ 0.02 — 0.08. Similar values have been obtained for l-D^I in Ref Q by 
considering the contributions of the charged higgs and top-pion to b — > sj. A similar exercise can be carried out 
with |-Dl6s|- Note that B s mixing probes the combination D* Lbs D m,bD* Rbs D lm, and so by either choosing R\ ~ or 
i?2 ~ we can satisfy the constraint on _B s mixing by choosing the appropriate mixing elements to be small. Note 
that in top color models we can have operators s(l — r y 5 )bd(l + 7 5 )d and s(l + j 5 )bd(l — r y 5 )d that can contribute to 
T — ► Bsd — > BX S after Fierz reordering. However these operators will be suppressed by form factor effects and also 
from mixing effects. We have checked that the contribution to T — > BX S from these operators are much suppressed 
relative to the contribution of the operators in L new . We will therefore not consider the the above operators in our 
analysis. 

Turning to R-parity violating susy we first collect the constraints on the relevant couplings. The upper limits of 
the L- violating couplings for the squark mass of 100 GeV are given by 



A', 



I A', 
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I A' 



1:5:5 



1 a; 



23j 



I A', 



3:5,, 



< 0.012, (k,j = 1,2,3;?: = 1,2), 
<0.16, (j = l,2), 

< 0.001, 

<0.T6, (j = l,2,3), 
<0.26, (j = l,2,3), 



(15) 
(16) 
(17) 
(18) 
(19) 



The first set of constraints in Eq. (15) come from the decay K — > -kvv with FCNC processes in the down quark sector 
fl3l| . The set of constraints in Eq. (16) and Eq. (18) are obtained from the semileptonic decays of i3-meson fl4|| . The 
constraint, on the coupling A' 133 in Eq. (17) is obtained from the Majorana mass that the coupling can generate for 
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the electron type neutrino [[L5|. The last set of limits in Eq. (19) are derived from the leptonic decay modes of the Z 
p6| . Assuming all the couplings to be positive we find the branching ratio for T — > BX S to be around 2 x for 
mp = lOOGeV. 

Turning next to J/ip — > DX U , we first make an estimate for this process in the standard model. Since the penguin 
c — > u transition is small in the standard model we neglect its contribution. As in the case for the T system , for a 
rough estimate, can write 

BR[J/i) -> D°X U ] ~ BR[J/tp -> D n }BR[D -> K-X]/BR[D° -> if "tt + ] 

We obtain BR[J/ip — ► D°7r°] from |3j and keeping in mind that £?i?[D° — » if contains contributions from states 
decaying to if" we obtain BR[J/ip -> - 1CT 10 . A similar exercise gives BR[J/tp -> D+ X u ] ~ 1CT 9 . 

Considering new physics effects we can constrain R\ and i?2 from c — > ul + l~ . We get an estimate of the constraint 
on c — > ite + e~ by adding up the exclusive modes 

Bi?[L> -> ue + e~] > Bi?[D -> (tt° + 77 + p + w)e+e"] 

From c — > ul + l~ one obtains 

|i?i, 2 |/A 2 < 3.7 x 10- 4 (l/GeF) 2 

We find the branching fraction for the process J /ip — » £>A U using the constraint from c — * ul + l~ can be (3 — 4) x 1CP 5 
In top color models taking R\ and R2 one at a time, one obtains 



2.1 x 10 3 

4 

or 

2.1 x 10 3 



ULcc\ 2 U Rt cUft t J 2 



2tt_ tt* |2 



1 \\URtc\ U Lcc U Lcu 

as the branching fraction for J/ip — > DX U . For mj between 100 — 200 GeV this rate can be between (0.1 — 2.0) x 10~ 5 
if all the mixing angles are *~ 1. It has been shown in Ref Ji"7| that our choices for and mj gives unacceptably 
large corrections to Z — > 66 from one loop contribution of the top pions. However in a strongly coupled theory higher 
loop terms can have significant contributions. Nonetheless if we change to ~ 100 GeV for better agreement with 
Z — > 66 data then the effect in J/tp — > DA„ is reduced by a factor of 16. As in the case of the T system we can satisfy 
the constraint from D mixing by choosing i?i ~ or R2 ~ 0. 

For R parity violating susy, contribution to J/ip — * DX U can only occur at loop level, with both the A' or A" 
contributing, through the box diagram and so is suppressed. 
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